Sequential magneto-tunneling in a vertical Quantum Dot tuned at the crossing to 

higher spin states 
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We have calculated the linear magnetoconductance across a vertical parabolic Quantum Dot with 
a magnetic field in the direction of the current. Gate voltage and magnetic field are tuned at the 
degeneracy point between the occupancies N = 2 and N = 3, close to the Singlet-Triplet transition 
for N — 2. We find that the conductance is enhanced prior to the transition by nearby crossings 
of the levels of the 3 particle dot. Immediately after it is depressed by roughly 1/3, as long as the 
total spin S of the 3 electron ground state doesn't change from S = 1/2 to S = 3/2, due to spin 
selection rule. At low temperature this dip is very sharp, but the peak is recovered by increasing 
the temperature. 

PACS number(s): 7I.30.+h, 72.I5.Gd, 75.30.Kz, 72.I5.Qm 



I. INTRODUCTION 

The conductance of a clean vertical Quantum Dot 
(QD) versus a gate voltage V g and a source-drain volt- 
age V s d allows the detailed study of the ground state 
(GS) and-first excited states of few electrons confined in 
the dottrel (see fig. |l|). The discrete levels due to the 
parabolic confining potential introduce a shell structure 
in the electron filling. In the case of a circular disk, ap- 
plication of a magnetic field B in the direction of the cur- 
rent (z-direction), orthogonal to the dot, favours states 
with increasing angular momenta, as well as higher total 
spin, also in the absence of Zeeman spin splitting which 
is believed to be small in these systems! The location of 
the crossings depends on the ratio between the Coulomb 
energy scale U and the level spacing of the confinement 
potential fluid- We consider here the case U ~ fluid- For 
larger ratios, transitions with B to higher spin states in 
dots with few electrons have been related to an "electron 
molecule" picture for the particle distribution within the 

dote. 

These features are strongly dependent on the long 
ranged interaction. Therefore, exactly soluble models 
like the Constant Interaction (CI) model or even the one 
with harmonic e-e interaction (HI) are unable to repro- 
duce thenxl. In fact, the HI model allows only for ground 
states which are singlet (S) for N even (N is the num- 
ber of electrons in the isolated dot) or doublet (D) for N 
odd. On the contrary, it was long ago recognized that a 
two-electron dot undergoes a Singlet-Triplet i£>-T) tran- 
sition whcrir-the magnetic field is large enoughcl. This was 
indeed seenH at B$t ~ 4 Tesla in dots whose confining en- 
ergy hud is about 5 meV. In fact, at difference with what 
happens in the 4 He atom (Bst ~ 10 5 Tesla), here the S-T 



energy separation at B = is only ~ 5 meV. Doublet- 
Quadruplet (D-Q) transitions for an N = 3 particle dot 
are correspondingly foundta (see below). Effects on the 
conductance are expected which are usually named "spin 
blockade" phenomenally: 

a) if the total spin of the ground state of the N + 1 
and N particle dot ( N+1 GS, N GS in the following) differ 
by more than 1/2, the dot is blocked, with the corre- 
sponding peak in the linear conductance missing at zero 
temperature. 

b) The reduction of the total spin in adding an extra 
electron to the dot has a lower transition rate as com- 
pared to spin increase. This can lead to negative differ- 
ential conductance. 

In this work we analyze in detail the behavior of the 
linear conductance of the dot when the gate voltage is 
such that the GS energies for N — 2 and N — 3 are 
degenerate (see fig. |l]b) and we show that, by increasing 
B beyond the S-T transition of the N = 2 particle dot, 
the conductance displays a sudden dip due to the spin 
selection rule b). Because the corresponding 3 GS has 
total spin 1/2, the conductance is depressed by roughly 
1/3. In fact, due to conservation of S z , it is impossible 
to add an electron with spin up (down) to the N = 2 
particle dot in the triplet state with S z = 1(— 1), to give 
total spin S = 1/2. These two (over six) possibilities 
that are missing are restored when the 3 GS changes to 
5 = 3/2. At low temperatures this dip is very sharp, but 
the peak is recovered in increasing the temperature. 

Three main ingredients are required to calculate the 
linear conductance, according to second order perturba- 
tion theory in the tunneling: 

i) the energies N E a of the isolated dot with e-e inter- 
actions (where N=2,3 and a labels the states of the dot). 
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studied point 

FIG. 1. (a) Sketch of the conduction band of the device in 
a linear transport experiment with a magnetic field B along 
the direction of the current, (b) Grey-scale conductance plot 
in the plane (V g ,V a d)- The circle encloses the conductance 



peak at the degeneracy point E 
studied in this work. 



(JV=3) 



(V g ,B)=Ei N = 2) (V a ,B) 



ii) the spectral amplitude Z N+1 '° for electron addition 
to the dot. 

Hi) the one-particle tunneling matrix T from the con- 
tacts to the dot. 

The first two quantities are provided by exact diago- 
nalization of the isolated dot, which uses a Lanczos code 
(up to 7 electrons), already employed previously in the 
analjesis of the magnetoconductance data in pillar struG- 
turestil. The third quantity was calculated analyticallylia. 
In Section II we study the transition rates in a dot from 
N = 2 to N — 3 by addition of an extra electron in 
an orbital of given quantum numbers. We analyze the 
various contributions to the tunneling, and the details of 
the spectral amplitude around B$t- We find that the 
spin selection rule mentioned above is responsible for a 
marked difference in the peak heights for transitions from 
2 T X -> 3 D 12 and 2 T X -> 3 Q 3 (see table I for the nota- 
tion). The linear magnetoconductance displays a sharp 
dip at low temperature at the S-T transition. In Sec- 
tion III we comment the results: similar peculiarities of 
the energy spectrum around Bst can be found for larger 
N (N > 3). They are remnant of the S-T transition in 
the N = 2 particle dot. However, in increasing N, the 
number of crossings is expected to increase and they oc- 
cur in a smaller energy range, so that their effect on the 
transition rate could be averaged out by temperature. 



II. LINEAR MAGNETOCONDUCTANCE AT 
THE S-T TRANSITION 

Transitions of an isolated dot to higher spin states 
with increasing magnetic field, follow from an interplay 
between the e-e interaction and the confining potential. 
They are easily monitored with exact diagonalization of 
the electron system. In fig. 2 the position in energy of 
the first levels of a two-dimensional(2D) QD for TV = 2 



and N = 3 is plotted versus magnetic field, expressed in 
Tilj c . Here oj c — eB/m*c is the cyclotron frequency and 
to* = 0.067m e is the effective electron mass in GaAs. 
States are labeled by the total spin S and by the compo- 
nents (along the direction of the field B) of the total spin 
S z and of the total angular momentum M (see Table I). 

The dot is laterally confined by a parabolic potential 
with level spacing fwJd = 4 meV. The e-e interaction 
is U x Zrf/|r — r'|, where Id — ^Jh/m*LUd is the length 
scale and the vectors are in the dot plane. The en- 
ergy scale of the Coulomb interaction U is taken equal 
to 3 meV. Spin splitting is not included in the picture 
for clarity. It is g* ^bB ~ 0.1 meV (/is is the Bohr 
magneton) at these fields if the bulk gyromagnetic ra- 
tio io^JGaAs (g* = —0.44) is assumed, but it could be 
lowerEj. A Lanczos algorithm is used on a basis which 
includes Slater determinants constructed with up to 28 
single particle orbitals that are solutions of the 2D har- 
monic confining potential. 

The S-T crossing in the 2 GS (fig. 2a,c) occurs at Uui c ~ 
7.5 meV. The next available levels are also shown and 
are located rather higher in energy. Correspondingly, 
the 3 GS (fig. 2b,d) is the doublet 3 D 1 at lower huj c 
and becomes the doublet 3 D 2 at hui c — 7.8 meV, ac- 
cording to the general rule that higher magnetic fields 
favor larger angular momenta. Next, there is a doublet- 
quadruplet (D-Q) transition at the magnetic field Bdq 
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FIG. 2. Energy levels of an isolated dot vs hu c . Zeeman 
spin splitting is not included, (a) and (c): dot with 2 elec- 
trons, (b) and (d): dot with 3 electrons. Fig.s (c) and (d) 
magnify the rectangle areas of (a) and (b) respectively. The 
level spacing of the confining potential is Tiujd = 4 meV, the 
strength of the Coulomb interaction is U = 3 meV. The GS 
levels 2 S, 2 T, 3 D 1 , 3 D 2 , 3 Q 3 are represented by solid lines (no- 
tation in Table I). Higher excited levels are represented by 
dashed lines. 
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FIG. 3. Transition rates as a function of uj/ujq where 
hid — E — 2 Eo is the energy difference with respect to the 
2 GS energy (wo = \J lo\ + u 2 /4). Zeeman spin splitting 
is not included. The step in Tiu> c is equal to 0.05 meV. 
Curves are shifted for clarity. A thermal broadening of 
T ~ 200mK has been added. The dot parameters are 
hujd = 4 meV, U = 3 meV. 



(huj c ~ 9 meV). Therefore, the sequence of transitions 
with increasing B, for addition of one electron at the low- 
est energy cost, is the following: 2 S — > 3 D 1 , 2 T ~ > 3 D 1 , 
2 T -> 3 D 2 , 2 T -> 3 Q 3 . The -> 3 £> 2 transition 
always occurs at a higher magnetic field than the S-T 
transition in the 2 GS: this has been numerically verified 
by varying the confinement of the dot. 

We calculate the transition rates for addition of one 
extra electron (nmcr) from the contacts, n and m are 
the principal quantum number and the angular momen- 
tum which label the one-particle energies e nm = hu>o(n + 
1) + mhuj c /2 (luq = y / uj 2 l + uj 2 /4). Here m has the same 
parity as n, taking the values —n, —n + 2, n — 2, n. 



The transition rate reads: 
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N, ) is the spectral 



weight amplitude. The dot is assumed to be in the N GS, 
\N0), prior to electron addition and a = (S,M) labels 
the final state with N + l electrons. 

The matrix T™^ appearing in eq. contains the sin- 
gle particle orbital overlaps. It describes the coupling 
of the dot to the contacts which have highly degener- 
ate Landau subbands and is factorized in the z and or- 
thogonal directions. It is also diagonal in the m in- 
dices, as the component of the angular momentum in 
the z— direction is conserved in the tunneling due to the 
assumed cylindrical symmetry. Peaks in the transition 
rate of eq. |l| are present each time the energy e matches 
Mat q 1 '" = N+1 E a — N Eq, provided conservation of S, S z 



and M are assured when adding the extra electron. We 
neglect Zeeman spin splitting: its effect will be discussed 
shortly in the conclusion. The transition rates are plot- 
ted in fig. ^ vs the dimensionless ratio oj/loo for increasing 
magnetic field. Here huj = E — 2 Eq is the excitation en- 
ergy with respect to the 2 GS. Curves (at intervals of 
Ahu) c = 0.05 meV) are shifted upwards for clarity with 
increasing lo c . An extra width has been given to the 
peaks (which are S— like in eq. Q), what is consistent with 
a thermal broadening of 200?™^ (see below). Each se- 
ries of peaks corresponds to the addition of one particle 
of definite angular momentum to the 2 GS . As an ex- 
ample, the peak from the singlet GS at cj/cjo ~ 2.3, for 
the lowest magnetic field, corresponds to the addition of 
one electron with m = 1, what changes the state of the 
dot from 2 5' to 3 D 1 . The initial and final states of the 
transition are labeled within the figure. The landscape 
changes drastically when the 2 GS becomes a triplet. 

In the plot of fig. || only the diagonal term n = n' 
has been included in the sums of eq. § This is expected 
to give the largest contribution within weak couplings. 

The spectral weight I^S^'™! is largely responsible for 
the heights of the peaks in fig. 0. 

This is seen in fig. ^a,b,c for the case when an elec- 
tron n = 2, m = 2 is added to the 2 T\ state (denoted 
symbolically as C\ 2a \ 2 Tl >), producing the final state 
with M = 3. There are four low lying states in the 
N=3 energy spectrum, labeled 3 Q 3 and 3 D 3 (fig. |]a). 
The lowest is the quadruplet state 3 Q 3 . The spectral 
weight for addition of an electron with spin down is dis- 
tributed among them (r.h.s. of fig. |Ja and fig. |i|c), while 
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FIG. 4. (a): The energy levels with M = 3 vs hu c are 
plotted on the right side. The spectral weights are depicted 
as horizontal lines to the right (+|) and to the left(+T) of the 
vertical line. (b),(c),(d): Spectral weig hts for the 2 T 1 -> 3 Q 3 
transition at 7lu;c=10.0 meV. huj = (E — 2 Eo) are the excita- 
tion energies. Dot parameters as in fig. 3. 
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in the case of a spin up particle the spectral weight is 
practically exhausted by the transition to the quadru- 
plet state with S z = 3/2 (l.h.s. of fig. |a and fig. |b). 
Up apd down arrows are to be interchanged in the case of 
2 Tl 1 ll3. Transition rates are proportional to the square 

of the Clebsch-Gordan coefficients \(S' S Z \SJL^ ± ^)cg| 2 
which accounts for the spin selection rulesEl 

These results are fully obvious in a picture in which 
Coulomb interaction is treated as a perturbation. 

Starting by filling single particle levels of the parabolic 
confinement potential, the N = 3 state (lowest in energy) 
with S z — 1/2 and M = 3 is obtained from the linear 
combination of three degenerate states with two electrons 
up and one down, distributed between n = 0, m = 0; n = 
1, m = 1; n = 2, m = 2. Coulomb interaction breaks the 
degeneracy, giving rise to the quadruplet state | 3 Q 3 / 2 ), 
which is the properly antisymmetrized combination of 
the previous three (with normalizing factor 1 / and 
two doublet states 3 -D 3 / 2 , which are almost unsplitted 
(see fig. ^a). The third doublet state at higher energy in 
fig. ||a is obtained by placing one electron in n = 3, m = 3 
and the two others in n = 0, m = 0. The transition 2 T\ 
— > 3 Q\j 2 by addition of one electron of spin J, has the 
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FIG. 5. Evolution of the height of the conductance peak as 
a function of 7kj c .(a) At a temperature of 200 mK, a dip is 
clearly visible in the range Tiuj c = 8-8.8 meV. (b) When the 
temperature increases, the dip leaves place to a peak which is 
due to the contributions of the 2 T 1 -> 3 D 1 and 2 T l -> 3 Q 3 
transitions. 



This spectral weight corresponds to the first peak of 
fig. |I]c. Similarly, the transition 2 T 1 — > 3 Qf/ 2 nas a 
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Comparing these noninteracting estimates with the ex- 
act results, the agreement is found to be close. 

To study the linear magnetoconductance we use a very 
general expression for the cwrent derived in the tunnel- 
ing Hamiltonian formalism^ If the coupling between 
the left (L) and right (R) barrier is proportionate, the 
linear conductance can be written in terms of the transi- 
tion rates of eq. [I]: 



9N,o(eV g ) 



dl 
dV 



N,0 



de t N Q (e) 



df 
de 



where / is the Fermi function. Just the diagonal terms 
are retained in eq. ^(what is found to provide the largest 
contribution by weak coupling, as above mentioned). 
Peaks are present when the gate voltage with respect to 
/1 is such that eV g = MatI) 1 ^ n the absence of crossings, 
only the GS (a = 0) contributes to the linear conduc- 
tance. Indeed, other a's would correspond to a nonequi- 
librium condition for the dot. Close to the S-T crossing, 
the conduction will be the sum of the two contributions 
gs and gx, weighted by the canonical equilibrium prob- 



ability of the dot being in the singlet or the triplet 2 GS 
(P s and P£ respectively, with ( = S z = ±1, 0)0: 



dl_ 

dV 



(eV 9 
v sd =o 



= Ps-9s+ Y, 
C=±i,o 



(2) 



In fig. [|, the maximum of the linear magnetoconduc- 
tance calculated using eq. ^| is shown at increasing mag- 
netic fields for two values of the temperature (fig. 3a: 
T = 200mK, fig. 3b: T = 600mK). At the lowest tem- 
perature it displays sharp features close to Bst and a 
marked dip between B$t and Bjjq. Conductance is at 
a maximum when the 2 GS is a singlet while the 3 GS is 
a doublet. The same happens at B > Boq when the 
2 GS is a triplet while the 3 GS is a quadruplet. In the 
latter case, there are six tunneling channels correspond- 
ing to 2 GS having S z — ±1,0 and addition of one extra 
electron with up or down spin. Each of them contributes 
1/6 to the total conductance. Between Bst and £?dq, 
however, the 3 GS is still a doublet, what inhibits two of 
the possible tunneling channels. 

Increasing temperature (fig. [5^3) , there is some overlap 
between the peaks due to the two doublet states 3 D 1 and 
3 D 2 and the state 3 Q 3 as well, which produces larger 
conduction in the intermediate region between Bst and 
Bdq and a new peak appears. This eventually washes 
out the dip. 
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III. CLOSING REMARKS 

The peak height in Coulomb blockade oscillations ver- 
sus gate voltage V g is governed by the N GS — ► N+1 GS 
transition rate. Using the S-T transition for the N = 2 
particle dot as an example, we have shown that, in the 
neighborhood of B$t, crossing of levels are to be ex- 
pected for the N + 1 particle dot as well. This feature 
seems to be general, as seen from the dot spectrum with 
N = 4 and N = 5. Increasing the ratio of the Coulomb 
energy U to the level spacing Tiujd due to the confinement 
potential, moves these crossings to lower magnetic fields. 

We have calculated the transition rates for addition 
of one electron to the dot with N = 2 in an orbital 
with quantum numbers nma. They are non-zero at the 
discrete energies e = /Ujvo While the single particle 
overlap T(e) plays a minor role, transition rates are deter- 
mined in weak coupling by the spectral amplitude for the 
isolated dot, which we obtain by exact diagonalization of 
a 2D system with parabolic confinement and Coulomb 
interaction between the electrons. 

Previous numerical work on the spectral weight of a 
few electron dot was mostly concerned with the .-total 
quantity w a = £ nm(T \Z^ a \ 2 (seS for B = andN for 
a = at various -B), in which the spectral weights are 
summed over all possible orbitals for the added electron. 
They showed that at definite energies its value could be 
strongly reduced with respect to the noninteracting case, 
due to correlations. However, in a vertical configuration 
for tunneling, with B in the direction of the current, tran- 
sition rates for the dot from the ^GS to the W+1 GS, or 
closely lying excited states at N + 1, require conservation 
of the total angular momentum M in the direction of the 
field, as well as of the total spin component S z . 

The quantity w a is then of little use, if conservation 
of the quantum numbers mentioned before has to be en- 
forced. We have shown that each of the contribution to 
w a , i.e. each of the spectral amplitudes, is strongly de- 
pendent on the quantum numbers which label the dot 
states and on the single particle orbital of the electron 
tunneling onto the dot. 

The spin selection rule we have been discussing in Sec- 
tion II arises from conservation of S z in a vertical ge- 
ometry. The absence of the quasiparticle (QP) peak for 
addition of an electron of a given spin shows that there 
are some tunneling processes that cannot contribute to 
the conduction. This produces a dip in the linear mag- 
netoconductance, in the range of magnetic field values 
for which the 2 GS — > 3 GS transition is between a triplet 
(5 = 1) and a doublet state (S = 1/2). 

As long as the Zeeman spin splitting is ignored, this 
mechanism is not spin sensitive, because with £ = 
0, ±1 are degenerate in energy. However, if Zeeman spin 
splitting becomes relevant, the 3 Q\/ 2 state could become 
lower in energy than the 2 D 2 j 2 and the range of magnetic 
field values between B$t and Buq, in which the dip in 
the conductance occurs, could disappear. 



Increased correlation, by distributing weight out of the 
QP peak to nearby energies, could reduce the effect. This 
is not the case, however, because we find that much larger 
U values do not change fig. |] qualitatively. Instead, they 
shift the magnetic field at which the crossings take place, 
and, of course, increase the energy differences between 
the levels. We do not find a sizeable spreading of the 
spectral weight to higher energies for U values as large 
as U = 15 meV. This could be partly due to the lim- 
ited number of orbitals used in the single particle basis, 
however. 
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level 




notation 


N=2, S= 


=0, M=0 (singlet) 


ago 


N=2, S= 


= 1, S K , M=l (triplet) 


2r,-,l 


N=3, S= 


=1/2, S x , M (doublet) 




N=3, S= 


=3/2, S z , M=3 (quadruplet) 


3 Qk 



TABLE I. Notations for the discrete levels of the QD used 
in the text. 
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